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Abstract 

We explore possible signals of CPT violation in neutrinos in the complete three-flavor framework. 
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describing non-standard interactions of neutrinos into the bounds on CPT violating quantities. 
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I. INTRODUCTION 



There have been theoretical suggestions that Lorentz invariance may not be an exact 
symmetry In such a case, even if the invariance is broken at a very high energy scale 
(say at the Planck scale in quantum gravity theories), the breaking is expected to leave its 
signature, however small, at laboratory energies. Such Lorentz violation may manifest itself 
as CPT violation. Indeed, in local field theories, CPT violation implies Lorentz violation 

In theories with spontaneous CPT violation [4J, the Lagrangian for a fermion to the lowest 
order in the high scale can be written as 

C = i^d^ij) - mijip - Ayfrf^ - B^zrf^ , (1) 

where and are real numbers. The terms containing A^ and B^ are clearly Lorentz 
violating, and give rise to an effective contribution to the neutrino Lagrangian that can be 
parametrized as 

C C V PTV = vlbf ^ 4 . (2) 

Here are four Hermitian 3x3 matrices corresponding to the four Dirac indices fi, wherein 
a, j3 are flavor indices. Then the effective Hamiltonian for ultra-relativistic neutrinos with 
definite momentum p is 

MM t „ 

H.^ + b, (3) 
where h = bo and M is the neutrino mass matrix in the CPT conserving limit. Following 



[i 



I], we choose to work in the preferred frame in which the CMBR is isotropic, where the 
rotational invariance implies no directional dependence for b. 

The same effective Hamiltonian can also be obtained by considering a modified dispersion 
relation for neutrinos, E 2 = F(p,m), in the presence of Lorentz violation. This dispersion 
relation may be written, using rotational invariance in the CMBR frame and demanding 
Lorentz invariance at low energy, as jlj] 

E 2 = m 2 +p 2 + E Pl fV\p\ + f^p 2 + ^\p\ 3 + ■ ■ ■ , (4) 

where /^'s are dimensionless quantities. The Planck energy Epi is introduced since it 
is the energy where Lorentz invariance is expected to be broken in quantum gravity. For 



ultra-relativistic neutrinos with fixed momentum p, the dispersion relation becomes 



2 

m 



E = P+^ + b + --- , (5) 

such that b = Epif^/2 is the leading CPT violating contribution. Generalizing this to 
three flavors leads to the same effective Hamiltonian as in ([3]). 

The possible origin of CPT violation in the neutrino sector has been studied in the 
context of extra dimensions jfJB], non-factorizable geometry Q], and non-local causal Lorentz 
invariant theories 9j. Bounds on the CPT violating parameters have been obtained in many 
different contexts. For example, the analyses of neutral meson mixings give |m^o — m-^o\ < 



1CT 18 m avg [nj], and \m B o — m-g-a>| < 1.6- 10~ 14 m avg ll|, whereas experiments on anomalous 
magnetic moment of muon put the bound on the anomalous frequency as — u)% | < 
10 -23 m M [3]. However, it is difficult to compare these bounds directly with the bounds 
obtained from the neutrino sector since we do not have an all-encompassing theory of CPT 
violation. 

The formalism to analyze CPT violating effects on neutrino oscillations has been proposed 
for the two-flavor case in The CPT violating contribution to the Hamiltonian would 
change the effective neutrino masses, which in turn would affect the neutrino oscillation 
wavelengths. Henceforth, for neutrinos we shall use p ~ E, so that M 2 /(2p) — > M. 2 /(2E). 
The typical frequency of neutrino oscillations is Am 2 /(2E), which can be as small as 10~ 22 
GeV in the atmospheric and long baseline experiments. Since the experiments measure the 
oscillation frequencies to an accuracy of ~ 10%, it may be naively estimated that neutrino 
experiments would constrain the CPT violating parameters to the order ~ 10~ 23 GeV. 

After the LSND result |13|] indicated three distinct neutrino mass squared differences 
when combined with the solar and atmospheric neutrino observations, it was proposed that 
the CPT violating effects may be large enough to make the neutrino and antineutrino spec- 
tra significantly different Q, [ijj] . However, this scenario was found not to be viable when 



combined with other neutrino experiments [16|, and the subsequent observation of oscilla- 
tions corresponding to AttIq in antineutrinos at KamLAND [17| ruled it out. If the LSND 
results are ignored in the light of the negative results of MiniBooNE 18] that explore the 
same parameter space, CPT violation is not required to explain any neutrino oscillation 
data. However, the current uncertainties in the measurements of Awq and Am 2 tm , which 



are ^ 9'.t and - 1 1't respectively allow the possibility of CPT violating effects in 
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neutrino oscillations, which may be observed or constrained at the future high precision 
neutrino oscillation experiments. 

CPT violation in neutrino oscillations would manifest itself in the observation P{y a — > 
vp) 7^ P(i>p — > Va). However, when neutrinos propagate through matter, the matter effects 
give rise to "fake" CP and CPT violation even if the vacuum Hamiltonian is CPT conserving. 



These fake effects need to be accounted for while searching for CPT violation. The u„ 



V,, 



channel was explored in this context in 



201 ]. where it was pointed out that CPT violating 



signals could become larger due to resonant effects. Using a two-flavor analysis, it was 
shown that a long baseline (L = 735 km) experiment with a typical neutrino factory setup 
can detect a difference of the eigenvalues of b upto 5b ~ 10~ 23 GeV. In 2l| it was shown 
that, using the atmospheric neutrino data at a 50 kt magnetized iron calorimeter with 1.2 
T magnetic field, 5b ~ 3 x 10~ 23 GeV should be clearly discernible in 8 years. The solar and 
KamLAND data gives the bound 5b < 1.6 x 10~ 21 GeV in (22]. Ref. 23] showed that for 
a hierarchical neutrino mass spectrum, the upper bound for the neutrino-antineutrino mass 
difference that can be achieved in a neutrino factory is \rri3 — rh^\ < 1.9 x 10~ 4 eV. Global 



two-flavor analysis of t 
5b < 5.0 x 10~ 23 GeV 



re full atmospheric data and long baseline K2K data puts the bound 



All the above analyses have been carried out in the two-flavor approximation. Moreover, 



it has been assumed in 



and 



2l| that the mixing angles as well as phases of the unitary 



matrices that diagonalize M and b in ([3]) are identical. The analysis of 22J takes the mixing 
angles to be identical and considers two specific values of the relative phase between the 
two unitary matrices. These assumptions have been made solely to simplify the analytic 
treatment, and do not have any physical motivation behind them. Ref. 24| analyzes the 
two-flavor case in its full generality, putting no extra condition on the mixing angles and 
the relative phase. However, a three flavor treatment is needed in order to obtain reliable 
results, since a two-flavor analysis cannot account for the CP violating effects that may 
interfere with the identification of CPT violation. The addition of the third (electron) flavor 
also compels one to take care of the matter effects when neutrinos pass through the Earth. 

In this article we consider the possible CPT violating effects that appear through ([3]), 
when three-neutrino oscillations are considered in their full generality. We treat the effect 
of the CPT violating term as a perturbation parametrized by a dimensionless auxiliary 
parameter e = 0.1 and express the differences of the eigenvalues of the CPT violating b 
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matrix, the reactor angle #13, and the ratio Am|/Am^ m as some power of e multiplied 
by 0(1) numbers, so that a systematic expansion in e can be carried out. The survival 
probabilities of and u e (and their antiparticles) can then be written down as a power series 
in e in a transparent form. This allows us to identify the combinations of CPT violating 
parameters that contribute to these probabilities to leading order in e. We compare the 
signals in the channels — > z/^, — > and v e — > u e , v e — > v e to estimate the extent 
to which these CPT violating combinations can be constrained or identified in future long 
baseline experiments. 

Bounds have been obtained on parameters describing the non-standard interactions (NSI) 

291 . 1301 ]. using both oscillation and non-oscillation 



of neutrinos with matter 



25 



26 



27 



28 



experiments. We explicitly show how to translate these bounds into the bounds on the 
CPT violating parameters in our formalism. This will allow us to compare and combine the 
bounds from these two approaches to restrict new physics in the neutrino sector. 

In the paper, Sec. [Til gives the parametrization of the CPT violating part of the effective 
Hamiltonian in flavor basis using the perturbative expansion scheme. Sec. II III and Sec. IIVI 
give the probability expressions and possible signatures in the long baseline experiments in /1 
and e channels respectively. In Sec]V]we summarize the current constraints on NSI param- 
eters and translate them to the bounds on CPT violating quantities. Sec. I VI I summerizes 
our results. 



II. PERTURBATIVE EXPANSION OF THE HAMILTONIAN WITH THE CPTV 
TERM 

In the three neutrino oscillation scheme, (z/ e , v T ) form the flavor basis and (v%, U2, ^3) 
form the mass basis, i. e. the basis in which MIVP / (2E ) is diagonal, and these are related by 

v « = [UoL^i , (6) 

where a G {e, u, r} and i £ {1,2, 3}. Let {v\, u\, v%) be the basis in which b is diagonal and 
let this basis be related to the flavor basis as 

= , (7) 
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where x G {1, 2, 3}. Both U m and U& are unitary matrices. 1 

When neutrinos pass through the matter, the electron neutrinos acquire an effective po- 
tential V e = \f2GpN e due to their charged current forward scattering interactions, compared 
to the other two flavors. Here Gp is the Fermi constant and N e is the number density of 
electrons. For anti-neutrinos, the sign of V e is reversed. The effective Hamiltonian in the 
flavor basis is 

diag(0, Amln Am^ 



U f « U - 



2E 



■ Vl + M b + diag(V e , 0, 0) 



where 



H 6 = U 6 -diag(0,& 21 ,& 3 i) -Uj 



(8) 



(9) 



Here b±, b 2 and 6 3 are the eigenvalues of b and bn = bi — bi for i — 2, 3. The net spectrum 
of neutrino mass eigenstates is given by the eigenvalues of Mf. We term the unitary matrix 
diagonalizing Hj as U/. 

A general N x N unitary matrix is parametrized by iV(iV — l)/2 independent real 
quantities (angles) and iV(iV + l)/2 independent imaginary quantities (phases). In case 
of the neutrino mixing matrix, (2N — 1) phases can be absorbed by redefining lepton and 
neutrino wave-functions. 2 Thus we can parametrize the 3x3 unitary matrix Uo by three 
angles 6*12, #23 and 6*13 and one phase 5 cp . We write Uo in the standard CKM parametrization 
as 



Uo = U23(^3,0)-U 1 3(^13,5 cp )-Ul2(^12,0) = U OT M({%};M 



(10) 



where \Jij(9ij, Sij) is the complex rotation matrix in the i—j plane, whose elements [Ujj] pg 
are defined as 

cos 9 p = q = i or p = q = j 
1 p = q 7^ i and p = q 7^ j 



[^(9,5)} 



pq 



sin 1 



p = i and q = j 



(11) 



— sin 9e tS p = j and q 







otherwise . 



In refs. 



20 



21 



22 



24 1, Uo and U& are 2x2 matrices. In addition, and 21 1 take these two matrices 
to be identical. 

2 For Majorana neutrinos, only N phases can be absorbed if M needs to be kept invariant. However, N — 1 
more phases are irrelevant when M only appears through the combination MM*. 
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Once we have redefined the phases of the lepton and neutrino wavefunctions to get Uo in 
the form ffTUl) . the basis of neutrino flavor eigenstates is completely defined. The matrix Ub 
then needs three angles 0ft 12 , 0623, 0&i3 an d six phases for a complete parametrization, given 
by 



U 6 ({0 Mi };{0«}; = diag(l,e^ 2 ,e 



i<j>b2 „i<t>b3 N 



u 



CKM 



diag(e iahl ,e iab2 ,e iab ' i ) 



(12) 



where a b x, a b 2, a b 3 are the Majorana phases and will not have any contribution to 
Mf through M b . Hence U/ may be written in term of a total of six mixing angles 
(On, 023, 013, 0M2, 0&23, 0fci3) and four phases (8 cp , 6 b , 4> b2 , 063)- 



21 



22 



24j arise from the 



Present limits on CPT violation in the neutrino sector [20l . 
limits on the neutrino oscillation wavelength, which in the two-flavor case gets modified as 
Am 2 tm / (2E) — > Am 2 tm / (2E) + 5b. The bound on 5b is therefore governed by the uncertainty 



in Am^ m , which is ~ 10%. Motivated by this, we assume that 



&2i,&3i < 0.1 x Am 2 atm /(2E ) « 0.13 x 10~ 2i GeVY^ 



(13) 



where Eq is the typical energy scale of the experiment. This may be parametrized by 
introducing two auxiliary quantities e = 0.1 and Se = 10~ 21 GeV 2 / E such that 



6 2 1 = ^2iSe , &31 = ^3iSe ■ 



(14) 



Clearly, fax, fox are numbers of 0(1) or smaller. The mixing angle 9x3 and the ratio 
AmyAmjj are small quantities, and may be expressed in terms of powers of e as 







13 = eXi3 



Am 2 21 /Am 2 31 = e 2 ( . 



The current bounds on the mixing angles and mass squared differences [19[ set 



Xi3 < 1. 



C-3.0. 



(15) 



(16) 



The sign of ( is positive (negative) for normal (inverted) mass ordering of neutrinos. 

Using the formal representation of 9x3, Arn^x and the CPT violating parameters in terms 
of powers of e as given in eq. (fT4j) and eq. (|15p . the Hamiltonian H/ can be expanded formally 
in powers of e as 



2E 



CO i ,2p(2) 



(17) 



where H^ ' 1 ' 2 '' are functions of all the mixing angles, phases, mass squared differences, and 
eigenvalues of b. All the elements of H^ ' 1,2 '* are of 0(1) or smaller, and Hy has nondegen- 
erate eigenvalues. The techniques of time independent nondegenerate perturbation theory 
can therefore be used to calculate the eigenvalues and eigenvectors of H/ upto the required 
order in e. These can be further used to calculate the neutrino flavor survival or conversion 
probabilities when neutrinos travel through matter with a constant density: 

Paf3 = P{v a -> vp) = 

where mf/(2E) are the eigenvalues of H/. This approximation is valid for neutrino propa- 
gation inside the Earth as long as neutrino trajectories do not pass through the core, and 
neutrino energy is not close to the #13 resonance energy in the Earth matter. 



— I- 



~2E 



(18) 



III. CPT VIOLATION IN P w AND SIGNATURES AT NOVA 

The survival probability of muon neutrinos of energy E, after traversing a distance L 
through the Earth is given as 



= 1 - sin 2 2fl 23 sin 2 A31 



+ efc , 1 A sin 2 2^3sin2A 31 - C 2 ^- sin40 23 sin 2 A31) + 0(e 2 ) , (19) 



^31 

where we define the dimensionless quantities A31 = Am 2 1 L/(4i?) and A = Am 2 1 L/(4i?o). 
The first two terms in ffl9|) are CPT conserving, and describe oscillations with frequency 
governed by Am 2 tm . The subleading contribution at 0(e) is CPT violating. The quantities 
C\ t 2 are given by 

d = cos 20 23 -5 2 sin 20 23 , 

C 2 = B x sin 20 23 + B 2 cos 20 23 , (20) 

where 

£?x = (H b22 - H 63 3)/(e^Eo) , B 2 = 2Re(H b23 )/(e^ ) . (21) 

The quantities B\ and B 2 depend only on b and U& as 

B x = [X cos 20 6 23 -Y sin 20 623 cos 5 b ] , (22) 
B 2 = — [X sin 20b 23 cos 50 + Y cos 20& 2 3 cos 50 cos 5b + Y sin 50 sin 5b] , (23) 

8 



wherein 



X = (3 21 cos 2 9 bl2 - An cos 2 9 bl3 - (3 21 sin 2 9 bl2 sin 2 9 m , (24) 
Y = (3 21 sm29 bl2 sm9 bl3 . (25) 

The phases <p b i only appear through the combination 5(f) = cf> b2 — <f> b3 . The corresponding 
quantity for antineutrinos can be obtained simply with the substitution e — > — e in f[T9l) . 
The terms involving matter effects as well as CP violation are suppressed due to 9 13 and 
AitLq, and appear only at 0(e 2 ) and 0(e 3 ) respectively. CPT violation can thus be cleanly 
extracted from the asymmetry 

Pw ~ P& = 2e ( CiA sin 2 29 23 sin 2A 31 - C 2 ^- sin 49 23 sin 2 A31) + 0{t 2 ) (26) 
V A31 / 

if it is indeed of the magnitude allowed by the current bounds. However, one needs to be 
away from the # 13 resonance, which for the Earth matter density occurs for E res « 5-10 
GeV, since the enhanced value of #13 makes the expansion in powers of e invalid. 

Eq. (jl9| demonstrates that, though the CPT violating parameter space consists of three 
angles, three phases and two eigenvalue differences (3 2 \ , P31, the effective CPT violating 
contribution to the probability P w is much simpler and depends only on two combinations 
of these parameters, B\ and B 2 , to leading order. A consequence of this result is that 
measurements in the muon channel can only put bounds on the two effective parameters Bi 
and B 2 , and not separately on the angles, phases or eigenvalue differences. 

Since 9 23 ~ 7r/4, the C% term in (f26j) dominates over the other. Accounting for the 
1/L 2 fall-off of the neutrino flux, the signal due to this term is optimized when L takes its 
minimum value that is able to satisfy sin2A3i w 1. This calls for a (relatively) low energy 
experiment with L/E ~ 240 km/GeV. The NOvA experiment 3jJ, |32| with its L = 812 km 
baseline and the NuMI beam energy E w 0.5-4.0 GeV satisfies these criteria, and hence is 
well suited to look for CPT violation. The energy range of NOvA is completely below the 
#13 resonance energy, so the contamination from CP violating # 13 contributions to — P^ 
is minimal. We take E for NOvA to be 1 GeV, so that S Eo = 10" 21 GeV. 

We demonstrate the validity and limitations of the analytic expression ( Tl9l) in the left 
panel of Fig. (TJ where is plotted as a function of energy for the NOvA baseline for the 
current best-fit values of Am^, Am 2 tm , 9% 2 , 9 23 and # i3 over the NOvA energy range. We 
choose normal mass ordering and Bi, B 2 with opposite signs, which is observed to be one of 
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FIG. 1: The left panel compares the analytical and numerical results for eight randomly chosen 
combinations of CPT violating parameters (thin lines) that correspond to B\ = B 2 = —0.3. The 
gray (shaded) band is the analytic expression plotted with an energy independent error of ±0.04. 
We choose normal hierarchy, Am| = 7.92 x 10~ 5 eV 2 , 9 U = 34.08°, Am 2 tm = 2.6 x 10 -3 eV 2 , 



42.13°, #i3 = 0.089 and 5 cp = 0. The right panel shows A^{E) for 4 years of running at 



NOvA with each fi + and fi~ , with an incident flux of 10 21 pot yr . The errors shown are only 
statistical. The central red (hashed) band shows the contribution in absence of CPT violation 
when the parameters are varied over their current 2a ranges. 

the worst case situations while comparing analytical results with numerical ones. We choose 
Bi — B'2 — —0.3 and take eight randomly chosen sets of CPT violating parameters that 
correspond to these values of B\ and B 2 . The plot shows that an energy independent error 
of ±0.04 in P w can account for the error due to neglecting terms of 0(e 2 ) or higher, in the 
whole energy regime of interest. 3 
We choose a typical NOvA setup 



31 



321 ]. with the NuMI beam directed towards a 0.5 



kt "near" detector placed 1 km away, and a 25 kt "far" detector at a distance of 812 km. 
The detector is assumed to be able to identify lepton charges. The neutrino propagation 
through the Earth is implemented using a 5-density model of the Earth, where the density 
of each layer has been taken to be the average of the densities encountered by the neutrinos 
along their path in that layer with the PREM profile 33j. We take care of the detector 



characteristics using the General Long Baseline Experiment Simulator (GLoBES) 34j, |35 | 



A part of this systematic error should be cancelled out when we concentrate on P^ — Pp.fi, however we 
choose to use a more conservative estimate of errors. 
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The cross-sections used are taken from 36|, [37|, and the simulation includes an energy 



resolution of cte = 10% y/E, an overall detection efficiency of 80% for all charged leptons, 
as well as additional energy dependent post-efficiencies that are taken care of bin-by-bin a 
la GLoBES. We assume perfect lepton charge identification, and neglect any error due to 
wrong sign leptons produced from the oscillations of the antiparticles. 
In the right panel of Fig. [TJ we plot the asymmetry 

_ Nj?(E) N^(E) 

M E ) = jynear^) ~ {E) ' ^ 

where Ne (Ng) is the number of l~ (£ + ) observed at the near or far detector. Here the 
events observed in the near detector act as a normalizing factor, and help in canceling out 
the systematic errors due to fluxes, cross sections and efficiencies in each energy bin. Note 
that modulo these factors, A^ is equivalent to (P w — P^p) multiplied by a geometric factor 
of (L near /Lf ar ) 2 . For plotting, we have considered a running time of 4 years with each of fi + 
and n~~, with an incident flux of 10 21 pot (protons on target) per year. 

The right panel of Fig. [1] illustrates salient features of the CPT violating contribution 
to Afi(E). The central band corresponds to possible signals in the absence of any CPT 
violating contributions, where we have varied #23, #13 and 5 cp over the currently allowed 2a 
ranges and have allowed for normal as well as inverted mass ordering. We fix Am|, Am 2 tm 
and 612 at their current best fit values, since variation with these parameters is not expected 
to be significant. For illustrating the signal in the presence of CPT violation, we choose 
5 cp = 0, \Bi t2 \ = 0.3, and fix 623 and #13 at their best fit values. The figure shows that A^ 
depends on both the magnitude and relative sign of Bi, B 2 and also on the mass ordering. 
It can be shown that the effect of changing sign of Bx,i is the same as changing the mass 
ordering, as expected from eq. (1261) when 623 ~ 7r/4. 

In order to estimate the possibility of identifying the CPT violating contributions from the 
experimental signals, with the current uncertainties in the standard three neutrino oscillation 
parameters, we display the confidence level contours in Fig. [2j In the left panel of Fig. [2, 
we have marginalized over all the standard neutrino oscillation parameters. It shows that 
B2 can be bounded from NOvA observations to the extent 

l-Sal < 0.1 (2a) . (28) 

The data are relatively insensitive to B\. This is expected from the analytic expression in 
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FIG. 2: Confidence level contours in the B\—Bi plane. The red (solid), green (dashed) and blue 
(dash-dotted) curves give the 3a, 2a and la contours respectively. An energy independent error of 
±0.04 on has been taken into account. Both the figures use 612 = 34.08°, Am 2 tm = 2.6 x 10~ 3 
eV 2 , 0i3 = 0.089 and B 1 = B 2 = as the input values. The left fi gure marginalizes over 023 , A?n 2 ^ m 
and 5 cp . The right figure uses 023 = 42.13°. 



eq. ( 1201) and (I26p : since 623 ~ tt/4, the terms containing B\ are highly suppressed. However, 
if #23 is known accurately and differs from ir/4, the sensitivity to B\ is restored. This becomes 
clear from the right panel of Fig. [21 where we have marginalized over all parameters except 
#23, keeping #23 fixed at a non-maximal value of 42.13°. In such a case, Bi can be constrained 
to 



|£i|<0.28 (2a) 



Accurate measurement of the deviation of #23 from its maximal value 



(29) 



is essential for the 



above bound on B\. The bound on B 2 , however, does not depend on the improvement in 
the measurement of any other quantity. A similar analysis performed for inverted hierarchy 
gives virtually identical results. 

The limits obtained in (j28p and (1291) are bounds on specific combinations of elements of 
H&. They imply 



H 622 - H b33 = elBilS^ < 10-^ GeV , 
2Re(H b23 ) = e\B 2 \S Eo < 1(T 23 GeV . 



(30) 
(31) 



In the limit 6^23 = v"/4,/3 2 i = 9 bl3 = S(p 
the quantity bounded in the two-flavor analysis 



0, the quantity B 2 in fact reduces to 5b, 



20 



2 II ]. Our three neutrino analysis thus 



12 



identifies the quantity that can be constrained, and also demonstrates that the constraints 
can be quantified in a clean manner at a low energy long baseline experiment like NOvA. 



IV. CPT VIOLATION IN P ee AND SIGNATURES AT A NEUTRINO FACTORY 

The survival probability for an electron neutrino travelling through a uniform matter 
density, in the presence of CPT violation, is given by 

P ee = l-4e 2 x 2 13 ( A A _ 3 \ ) sin 2 (A e -A 31 ) 



31 



- 4e 2 xi3 [Re {S<*>)] Aai *° sin 2 (A e - A 3a ) 



2n _ x A 2 2A 2 - 2A e A 31 + A^ 



e Pi cos 5, 



cp 4A 2 (A e - A 31 , 



) 2 



A 2 , A 31 (-2A e + A 2 
'4A 2 (A e - A 31 ) 



+ e (P 2 cos 26*23 + P3 sin 26> 23 ) cos 5 cp — ^ — — — - 2 — 

1 + 0(e 3 ) , (32) 



+ e 2 



A 2 2A 2 

P 4 -n7Cos2A e - \S\ 2 — £-^cos (2A e -2A 31 ) 

2A 2 1 1 A e - A 31 2 v 4 ' 



where A e = V e L/2 and recall that 9\ 3 = exi3- The CPT violating quantities appearing in 
(J3"2"|) can be expressed in terms of two complex quantities Q = Hm. 3 and R = H&i 2 : 

Q = Q 1 +iQ 2 = -icos^ 13 e^ b3 |/? 2 isin2^ 12 sin^ 23 

-2 (/3 31 - (3 21 sin 2 6 m ) sin # bl3 cos e m e~ i5b J , (33) 

R = R 1 + iR 2 = l - cos e bl3 e-^ |/? 21 sin 2fl bl2 cos # fe23 

+2 (fti - 02i sin 2 612 ) sin 6 bl3 sin 623 e~ Wi J , (34) 

where Qi and Ri are real numbers. In terms of Q and P, the CPT violating parameters S 
and -Pi, 2,3,4 in (1521 may be written as 

5 = Q cos8 23 + Rsin9 23 , (35) 
Pi = |Q| 2 + |P| 2 , P 2 = \Q\ 2 - \R\ 2 , P 3 = 2Re(QR*) , P 4 = Pi - |S| 2 . (36) 



In eq. (1321) . the first two terms are the expression for P ee with a non-zero #i 3 when CPT 
is conserved, while all the other terms are CPT violating contributions. There are no 0(e) 
terms. Both the 9i 3 correction as well as the CPT violating contributions appear at 0(e 2 ). 

13 



0.94 
0.92 
0.9 
o 0.88 

1. 

0.86 
0.84 
0.82 
0. 



15 



20 



Q 1 =Q 2 = R 1 =R 2 = 1.0 

25 30 35 40 
E (GeV) 



45 



50 



300 
200 



» 100 



X 
^-100 
-200 
-300 



5 CD = ; NH 
p IH 
8 Pn = 7t/2; NH 
IH 

5™ = it ; NH 
IH 



-a -- E- S 5" -> ■ 



-J- 
15" 



20 



25 



30 35 
E (GeV) 



40 



45 



50 



FIG. 3: In the left panel, the black (solid) line is the analytic expression, and the gray (shaded) 
band corresponds to the analytic value with an error ±0.015. The lines in six different colors 
(symbols) are for six random sets of CPT violating parameters with Q\, Q2, R\, R2 fixed at 1.0. 
We choose normal ordering, and the values of Attt-q, #12, Am^ tm , #13 and 5 cp the same as 
that in Fig. [TJ The neutrinos traverse through the Earth for L = 3000 km before being detected. 
In the right panel, the central band shows the contribution in absence of CPT violation when 
the parameters are varied over their current 2a ranges. The counts are for 2 years of running of 
neutrino factory with each of e + and e~ . The errors are only statistical. 

These also include terms that get contributions from both #13 and CPT violating parameters. 
The probability P ee depends only on two complex combinations Q and R of CPT violating 
parameters, so this channel can put bounds only on these two parameters. 

Note that the coefficients of 0(e 2 ) in ( 1321) contain terms proportional to (A3i/A e ) 2 , which 
should be small for the e-expansion to be under control. Therefore, the expression fl32|) is 
valid only when A e > A31, which happens at energies above the #13 resonance energy. In 
order to get significant effects at large energies, one also needs long baselines. Both these 
conditions would be satisfied at a neutrino factory with an energy range 10-50 GeV and a 
baseline ~ 3000 km. If we restrict ourselves to energies well above the #13 resonance energy 
~ 5-10 GeV, even the CPT conserving 9f 3 contribution is suppressed, so that the CPT 
violating contribution can be more cleanly identified. For the neutrino factory, we can set 
the typical energy E = 10 GeV, so that S Eo = 10~ 22 GeV. 

In Fig.[3]we demonstrate the validity and limitations of the analytic probability expression 
f l32|) . where we choose the mixing parameters Am|, Am 2 tm , 6 ) 12 , #23, #13 to have their best-fit 
values 19|. We choose normal mass hierarchy, and fix S cp = 0. This is observed to be one 
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of the worst case situations while comparing the analytical expressions with the numerical 
ones. For the CPT violating part we choose Q\ = Q2 = R\ = R2 = 1.0, and six random 
choices of the elements of b that map to these values of QiS and RiS. It is seen that an 
energy independent error of ±0.015 on P ee can account for the error due to neglecting higher 
order terms in e over the whole energy range of interest. 

To demonstrate the capability of a typical neutrino factory setup for identifying the CPT 
violating contributions, we define the asymmetry 

A(F) - Nl-(E) N f :\E) 



in a typical neutrino factory setup [39] with a 50 GeV muon beam directed to a 0.5 kt "near" 
detector 1 km away, and a 50 kt "far" detector 3000 km away. The detectors are assumed to 
be capable of identifying lepton charges. The number of useful muons in the storage ring is 
taken to be 1.066 x 10 21 , which corresponds to approximately two years of running with pT 
and /i + each at the neutrino factory, using the NuFact-II parameters in {^(J . The simulation 
includes an energy resolution of (Je/E = 15%, and an overall detection efficiency of 75% 
for all charged leptons. Earth matter effects, interaction cross-sections and post-efficiencies 
are taken care of in the same way as was done in the case of NOvA. We assume perfect 
lepton charge identification, and neglect any error due to wrong sign leptons produced from 
the oscillations of the antiparticles. GLoBES is used to get the energy variation of the 
asymmetry A e (E) as shown in the right panel of Fig. [31 The figure indicates that it will be 
possible to discern the CPT violating contributions from the background of CPT conserving 
contributions if, for example, Qi = 0.8 and Ri = 0.5. Note that A e is approximately 
equivalent to (P ee — Pgg) multiplied by a geometric factor of (L near /Lf ar ) 2 . 

The magnitude of A e {E) depends on Qi, Q 2 , Ri, R2, 8 cp and mass ordering. To estimate 
the possibility of identifying any CPT violating signal in spite of our current lack of knowl- 
edge about the standard oscillation parameters in the CPT conserving case, we display the 



confidence leve 



contours in Fig. HI We have chosen the best fit values of Am^, #12, Am 2 tm , 
623 and #13 [19] as the input values. Since we do not have any information about 5 cp , we 
choose the input value of 5 cp in the range that is observed to give the most conservative 
bound on \Q\ and \R\. From the left panel of Fig. HI for the normal mass ordering the 
bounds obtained are 

|Q| 2 <1.1, |i?| 2 <1.35 (2a), (38) 
15 



while the right panel with inverted mass ordering gives 

|Q| 2 <1.2, |i?| 2 <1.4 (2a). (39) 

It is observed that if the actual value of 823 is smaller, the \Q\ 2 bound decreases and the 
bound on |i?| 2 becomes larger. The reverse is true when 823 value is higher than the current 
best-fit value. This is true for both the mass orderings. 




0.4 0.8 1.2 1.6 2 0.4 0.8 1.2 1.6 2 



|R| 2 |R| 2 

FIG. 4: Confidence level contours on |Q| 2 -|i?| 2 plane. The red (solid), green (dashed) and blue 
(dash-dotted) curves give the 3a, 2a and la contours respectively. An energy independent error 
of ±0.015 on P ee has been taken into account. We use the same AmL 612, |Am 2 tm |, 823 and #13 
input values as in Fig. [TJ The additional input values are 5 cp = ir/3 and \Q\ 2 = \R\ 2 = 0. All the 
parameters other than \Q\ 2 and |i?| 2 are marginalized over in the analysis. 

The bounds on \Q\ 2 and |i?| 2 translate to 

|H 613 | = e\Q\S Eo < 10- 23 GeV , (40) 
|H 612 | = e\R\S Eo < 10- 23 GeV. (41) 

The reach of A e for the CPT violating observables is thus similar to that of as obtained 
in Sec. IIII1 However, note that the actual combinations of elements of Hf, constrained by 
the muon and electron channels are quite different. 
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V. CONSTRAINTS FROM BOUNDS ON NON-STANDARD INTERACTIONS 



In the presence of NSI of neutrinos with matter, the effective Hamiltonian in the three- 
flavor basis becomes 



diae(0, Am,, , Ami i 
H ; « U • Sl 2E 21 ' ^ • Uj + V e e NSI + diag(K, 0, 0) , (42) 



where ejvs; is a 3 x 3 matrix 



e NSI 



^ee ^efi ^er 



e e^ e A i A t e M r 



(43) 



^ ^ ^ 

y t er t fir tt J 

that parametrizes the NSI interactions. The factor of V e multiplying €nsi represents that the 
net NSI strength depends on the density of matter. The Hamiltonian for the antineutrinos 
will be obtained just by V e — > — V e and t a p — > e*^. 

Since CPT violation necessarily implies NSI, the bounds on the NSI violating parameters 
e a /3 would restrict CPT violating parameters as well. In order to see the exact correspon- 
dence, note that the oscillation experiments are sensitive to only differences in the eigenval- 
ues of the hamiltonian, and not to the absolute eigenvalues. Therefore, the part of the NSI 
relevant for oscillation experiments is only 

IRnsi = £nsi - e ee I , (44) 

where I is the identity matrix. Then the comparison of eqs. (jHJ) and ( 1421) implies that the 
mapping 

M b M NSI (45) 

would allow us to translate the results from one parametrization to the other. Note that 
there is a difference between the two sources of nonstandard physics under consideration. 
Whereas Mnsi is proportional to the matter density, is independent of it. However, 
as long as the matter density relevant for the experiments restricting Ha^/ is known and 
is almost a constant, the CPT violating contributions may be mimicked by the NSI ones. 
Therefore, the bounds on e a p from the NSI analysis can be translated to the bounds on the 
elements of H& in the CPT parametrization. 

Note that the bounds obtained from the CPT analysis cannot be applied to the NSI 
bounds, since there can be sources of NSI that are CPT conserving. If an experiment is 
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sensitive to the variations of matter density along the neutrino path, it will be able to 
separate the NSI contributions from the CPT violating ones. 

A two flavor analysis of the atmospheric neutrino data combined with the MACRO data 
25( and K2K data [26J] yields 



0-05 < e MT < 0.04 (99% C.L.) => |H b23 | < lO" 23 GeV 



(46) 



where we have assumed an average density of 4.5 g/cc inside the Earth. This bound is 
comparable to what would be obtained using long baseline experiments as described in 
Sec. ED 

The neutrino scattering experiments CHARM and NuTeV mainly constrain the NSI 



couplings of z/ M , and give 



27 



41 



42 1 



M £ 10" 



H 623 | < 1.1 x 10" 25 GeV 



25 



(47) 



where we take the average Earth matter density to be 2.7 g/cc. This constraint is extremely 
strong, and would imply \R\ rs 0, thus simplifying the analysis of Sec. IIVL These experiments 
also bound 



Iwl £ !0" 



(48) 



which by itelf does not put any constraints on the CPT violating parameters since only the 
differences between the diagonal elements of new physics hamiltonian is relevant. 
Using the bounds on e Mj g stated above, 
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29 



30j analyzed the possibility of constraing 



e ee , e eT and e TT in MINOS experiment assuming e efJi 
two-neutrino analysis leads to 



K\ < 2.9 



C.L.) 



-fJ,T 



0. This effectively 



H 6 i 3 | < 3.2 x 10-" GeV 



(49) 



which will be improved significantly at the neutrino factory with the v e — ► v e channel, as 
described in Sec. HVl 

The 99% C. L. bounds on the diagonal NSI elements, given the initial assumption of 



e MM = 0; translate as 



29] 



- 0.4 < e TT < 4.5 -0.5 x 10~ 22 GeV < M b22 - H b33 < 5.0 x 10" 22 GeV (50) 
-1.0 < e ee < 0.9 |H b22 | < 10~ 22 GeV . (51) 
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Here, we have taken an average matter density of 2.7 g/cc, which is relevant for the MINOS 
baseline of 732 km. As seen in Sec. IHfl NOvA will be able to constrain H&22 — H&33 to a much 
better accuracy. The channels we have considered are rather insensitive to the absolute value 
of |H 622 |. 

Current and future long baseline experiments like OPERA and T2KK are expected to 
improve the bounds on NSI parameters 



43 



44|, and hence indirectly, those on the CPT 



violating parameters. Data from a future galactic supernova will also contribute to con- 
straints on NSI parameters 45|, but converting them to bounds on CPT violation will 
not be straightforward since the situation cannot be approximated with a constant matter 
density. 



VI. SUMMARY AND CONCLUSIONS 



We have calculated possible CPT violating contributions to neutrino masses and mix- 
ings in the complete three-flavor analysis. Parametrizing the leading CPT violating effects 
by a Hermitian matrix b that adds to the effective neutrino Hamiltonian, we have de- 
veloped a framework based on the perturbative expansion in a small auxiliary parameter 
e = 0.1. It involves expanding the elements of (the matrix b in the flavor basis), the ratio 
Amg/Am^ m , and #13 as powers of e multiplied by 0(1) numbers. This allows us to treat 
the CPT violating b contributions in all generality, while keeping the analytical expressions 
simple and transparent. Though the complete parametrization of b involves three eigen- 
values, three mixing angles and six phases, we show that only certain combinations appear 
in the survival probabilities of muon and electron neutrinos, so that the analysis needs to 
concentrate only on limiting those combinations. 

The survival probabilities of and to 0(e) involve only two combinations of elements 
of Hi, viz. the real parameters B\ oc H;,22 — and B 2 oc Re(H;,23). Formally, the 
CPT violating contribution due to these terms is of a higher order than the CP violating 
contribution in the CPT conserving limit. The contribution due to B\ vanishes when #23 is 
maximal, so that a deviation of 6*23 needs to be established in order to put any bounds on 
this parameter. The other quantity B2 may be constrained to be \B%\ < 0.1 with 4 years of 
running with and each at NOvA with an incident flux of 10 21 pot yr _1 at 2a. This 
would correspond to bounds on H&22 — H 33 and Re(H;,23) of the order 10~ 23 GeV. Note that 
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though the constraints that we have obtained are of the same order as those obtained in 



earlier studies 
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211 ]. our analysis identifies the exact combination of elements of H&, and 



hence b, that these bounds apply to. 

The CPT violating contribution to the survival probability of v e and i> e appears only 
at 0(e 2 ), and hence is expected to be more difficult to extract. We isolate two different 
combinations of elements of b, viz. the complex parameters Q oc H&13 and R oc H&12, that 
govern this contribution and numerically analyze the feasibility of extracting them. We 
demonstrate that for \Q\, \R\ > 1.2, it may be possible to ascertain the presence of CPT 
violation at 2a at a neutrino factory with a detector at L = 3000 km that can distinguish 
u e from i> e , within 4 years. This corresponds to bounds on Hbi 2 and of the order 10~ 23 
GeV. Note that the exact combinations of elements of that are constrained by the muon 
and electron channels are quite different. 

The CPT violating observables and A e in this paper are the same as those considered 



m 



46| for disentangling the signals of sterile neutrinos. The energy dependence of the 



signatures of CPT violation and sterile neutrinos, however, is different and these two new 
physics signatures may be disentangled with a combined analysis. 

The constraints obtained on the NSI parameters through oscillation and non-oscillation 
experiments can be translated to bounds on elements of H&. We find that the bound on | H512 1 
implied by the NSI constraints is much stronger than the expected reach of even neutrino 
factories, whereas the bound on IH&23I is comparable to the one expected at NOvA. On the 
other hand, | H513 1 and the difference H&22 — will be much better constrained by the 
long baseline experiments. NSI analyses give a constraint on the absolute value of H(, 2 2, to 
which the channels we have considered are rather insensitive. 

In this paper, we have confined ourselves to low energies (E < 5 GeV) for the muon 
channel and high energies (E > 15 GeV) for the electron channel. This allowed us to cleanly 
isolate certain combinations of elements of H&, viz. two real quantities H&22 — H^, Re(H;,23) 
through the muons and two complex quantities Hbi2,Hbi3 through the electrons. A more 
exhaustive analysis that uses the complete energy range and the long baseline as well as the 
atmospheric neutrino data may lead to constraints on other combinations of elements of Hft. 
However, it is not clear if it can be achieved through a clean analytic treatment. 
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